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Contribution to the Theory of Attraction when the Force Varies 

as any Power of the Distance. 

By Major R A. MacMahon, F.E.S., and H. B. C. Darling, 
(Eeceived April 8, 1918, — Revised December, 1918.) 

Introduction. 

The object of this paper is to make a contribution to the theory of 
attraction when the force is proportional to any given power of the distance. 
In the case of the law of nature, it is well known that when the attracting 
mass is a hollow shell of uniform density, whose exterior and interior 
bounding surfaces are both surfaces of revolution about the axis of z, the 
potential at any point exterior to the exterior surface can be simply obtained 
from the potential at any point on the axis exterior to the exterior bounding 
surface, and a similar theory connects the potential at any point within the 
inner bounding surface with that at any point on the axis which is also 
within the inner surface. In the expression for the point % on the z axis we 
have merely to write r for z, and intercalate the zonal harmonics P 0) Pi, Pa, . . . 
in the well-known manner. 

We consider the analogous theory when the law is that of the nth power 
of the distance instead of that of the inverse square. 

In a notation which we find the most convenient for the expression of our 
results we consider the expansion 

so that when the second suffix of P is zero the coefficients are zonal 

harmonics 

o, o> ri, o, X2, o 

instead of Po> Pi, P2, 

in the ordinary notation. 

These double-suffixed coefficients have been considered by Heine and other 
investigators,* but not in the present notation or with the present object in 
view. 

In the result we find that the expressions to be intercalated are linear 
functions of these double-suffixed coefficients, the coefficients which connect 
them being derived from potential functions in regard to a point on the z 
axis for different values of s. 

* Heine, ' Kugelfunctionen/ Band I, p. 297, 
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The differential equation satisfied by V m , g is 

(1— fi 2 )d f 2 'P m)S ^(2s— 2)^P m ,, + w(m — 28+l)J?. m , t =s 0. 

Integrating this equation with respect to p in the ease m = 2s, we obtain 

(1 — fjL 2 )d tl J y 28>s-\-%SfiT?2s, s = constant. (1) 

When s is an integer, the left-hand side of (1) vanishes for > == ; thus the 
constant in (1) is zero. When s is fractional, we may suppose P & , ,■ subject 
to the condition that the constant in (1) is zero. Hence, from (1), we obtain 

P*„ = C(1 -/*»)■, 

where C is a constant with respect to jjl. Considering the expansion of 

(1— 2fJLOL + 0L 2 y 



\S"-i 



when fi = 0, we see that C = 1, ^'^;"^ 2g X \ and therefore 

p _ 1.3.5 (2s— 1)^ 3 y /m 

NOW put Pw+2s, » — Pg*, * T m , s 

and substitute in the differential equation satisfied by P m +2*>*. We find on 
reduction 

(l-/x 2 )3/T w ,,-(2s+2)^T w , g f m(m+2s+l)T*„ - 0, 

which is seen to be the same as the equation satisfied by P, n , ~ s . Hence 
since T m s and P w , _, are of the same degree in jjl they can only differ by a 
numerical factor. 

where A is a numerical factor to be determined. 

Putting /jl = in this equality we find that when m is uneven both sides 
vanish, but when m is even 

( 2s-l)(2s~3) (-m + 1) 

2.4 (m + 2s) 

_ A 1.3...(2s-l) , v^ (2s± l)(2g + 3)......(2* + m-l) . 

2. 4.. ,2s " 2.4 m ' 

and on reduction 

1 = A /m+2 S \ 
\ m j 

Thence I j P OT +2*,* = P2*, *P?»» -•> (3) 

and therefore by (2) we see that if n > 2s— 1, then 

P n , , == 0, for /* == 1. 
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A knowledge of this property is necessary in what follows. The 'relation 
(3) may be deduced from results already published which have been other- 
wise obtained by different authors. 

1. In trie application of zonal harmonics to physical questions it is usual to 
deal in the first instance with surfaces and solids of revolution symmetrical 
about an axis. In the theory of attraction according to the law of nature we 
have in such a case an easy way of finding the potential at any point which 
is not inside the attracting material when once the potential at any point 
upon the axis of symmetry is known. As usually stated :— 

Suppose that an attracting mass, M, of uniform density, is a hollow shell 

whose exterior and interior bounding surfaces are both surfaces of revolution, 

their common axis being the axis of z ; let the origin be taken within the 

interior bounding surface, and let E be the distance from any point of the 

mass to a point upon the axis of z at a distance z from the origin and within 

the inner surface; then if we write the potential at the point upon the 

axis of z 

^K^dM == M2a», 3 n 

the potential at another point, not upon the z axis but within the inner 
surface whose polar co-ordinates are r,9 may be expressed by 

M%a n> o P n (cos 6) r n . 

And there is a similar theorem with regard to points which are without the 

outer surface. 

The problem before us is the determination of the analogous result when 

the potential is of the form 

XR 2s ~ l dM 

where s may be any real numerical magnitude. 

The result is found to involve as coefficient of r n a linear function of the 
coefficients P m , s . 

2. The effect of Laplace's operator 

V 2 

upon the product P myS r n has been given by other investigators in a different 
notation. 

The formula in the present notation is 

V a Pm,,r»= {(w-m + 2s)(^ + w-2« + l)P TO , J r + 2c9(2s--l)P mil -i}r»-- a 
leading to the particular case 

V 2 P»,,r* = 2s(2s-l')P n - 3 , f - 1 r»- a . 

3. We first require the following Lemma :— 
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:f The value of the integral 

'2tt 

Pn> s {cos cos 0' -f sin sin 0' cos (<£ — </>') } ^^ 

, 

may be always expressed in the form 

GoP n , o (ft) + CiP n , l (/*) 4- C 2 P n , 2 (/<0 + . . . 

where /^ = cos ; c , ci, e 2 , ..., are independent of /a] and the series is carried 
to I (n + 1) or | (% + 2) terms/' 
For, since all uneven powers of 

sin sin 0' cos (<£ — <£') 

will disappear on integration, we see that the integral will be an expression 
of the nth degree in /j,, the successive powers of fi being n, n~~2, n — 4, ...... 

Moreover, this is also the case with the series 

c V n} o (/a) + o x B n , i (fi) -jr . . . 

and the unknown quantities c , ci, c a , ... are just sufficient to determine the 
equality without ambiguity. Further, the functions P which occur in the 
series will all be of the form V mt where n>2t — l i and will therefore all 
vanish when p = 1, except for t = 0. 

Hence when fi = 1 the series vanishes with the exception of the first term, 
which is 

c P nj o<X) = Co. 

We now write 

r2w 

r n P«, « {cos cos 0' + sin0sin0'cos(<£ — §'))cl$ f 
Jo 

= ^ {Co'P», O) + CiP„, i (/a) + C 2 P WJ 2 (/&)+... } (1) 

and putting //, = 1, we have 

'2tt 

Pn, s (COS / )d(j) f = C 



or Co = 27rP n , ,(///) where /a' = cos 0', 

Now applying Laplace's operator 

'27T 

2 s (2s- 1) r*"* 2 P w - 2 , s _! {cos cos 0' + sin sin 0' cos (<£ - <£')} ety' 

Jo 

= r"~ 2 {l . 2 cxP w .j,o(/a) + 3 . 4caP B - a , 1 (^) + 5 . 6 c 3 P„- 2 , 2 (//,) + ...} 
and herein putting /a = 1, we find that 

C! =( 9 / . 27T-P„-2,.-l(/Lt / ). 



Ot/4: IVLsblOi Jr. J\j,t i^lM3J3xcL 003.1 cbOQ. ifULx. JlX,« Jl5. v*/» XJaixXkJjLQ.* 



Again applying Laplace's operator we find 



and generally 



c 2 



y p 



'2 s" 



j | . 27Tx n — 4, g»~%\fM ) 



/'2s\ 

2.W 



. 27rP n -2 



5jS>5 «~"^ 



(/A 



'2 a 



where we must remember that l%° ) is an abbreviated notation for an 

expression in which s may be any real number. 

Hence putting 

cos cos 6 f ~f sin sin 0' cos (</> --<£') = cos © 
we see that 



J2ir 



P», s (COS ©) rf0 ; 



2s N 



= 2 7r {P w , , (^') P M , o (/*) + ( 2 ) P*-* *-! (/*') -P»» 1 (/*)+•"•• • 

+ ( 2j0 ) Fn-^.-*^') *»*(/*) +.•-}, (2) 

an expression in which p, p may be interchanged. 

4. Now in the case of any solid or surface of revolution, the potential, 
which we will write V„ may be expressed respectively in the forms 



f 



r2n 



/(?•', 6') (r 2 - 2rr' cos © + r /2 ) s ^ d$ dff dr' 



(3) 



and 



'2w 







f(r\ ff) (f-2rr' cos © + /*)•-* ety' d0', 



where cos © has the meaning assigned above ; r, 0, $ are the co-ordinates of 
the fixed point for which Y s is taken ; and f(r',0') is some function of the 
co-ordinates r',0' of a point (r'0'0') of the solid or surface respectively. 

Our remarks will apply mutatis mutandis to both solid and surface, so that 
for brevity we restrict ourselves to the triple integral above set forth. 

First suppose r to be less than r' throughout. 

It may be given the form 



r *f2ir 
J J J » 



(4) 



For a point upon the z axis distant z from the origin, so that = 0, we may 
suppose this expression equal to 



JXJl. ^fWfi f gZ , 



so that 



Ma, 



m 



r'*-«- 1 /(r\ 6') P., , 0*') «!f d&' dr', 
= 2wJJr'*-»- 1 /(r', 0') P n , ,(/*') d0' dr' ; 
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whence also 

and in general 

Ma*-*,-, = 2tt jjr' 2 *~»~ \/(/, 0yP n ^,,-,(/)^'tfr'. 

Hence by (2) we see that the expression (4) for the potential may be 
written 

M 2 -| a n> t P n , o (/u,) -f a n ~2> «- i ( 2 * ) &n> i <» + 



71 ss QO 



'2s 



or MS 2( 2 J^-^t-j»Pn, ?(/*)*"; ( 5 ) 

in which the upper limit of p is reached when n — 2p in either unity or 
zero ; so far as the series does not terminate automatically before this point 

is reached by the vanishing of ( ~ s 

5. Next suppose r>r'. 
Then (3) may be written 

f f [*/(/> V) 2 r'»r~<»- 2 * +1 >P n , s (cos 6) d<f>' dff dr'. (6) 

For a point upon the z axis, where 6 = 0, we suppose this expression 
equal to 

and then we have 

mab " = fffr r "' /(r '' ^)p».»(m')#'^'**' 

= 2 W ||r'V(r',0')P»,.(/)^'^' 

leading to the general relation 

MA„, P = 2tt Hr'Hfir',?)^, „((*') d&dr' 

and therefore by (2) 

SMA B , , (l s \ P B _2p, .-, (?) =* 2,r JJ /»/('•', *') Pn, . (cos @) «W d/, 

and the potential becomes 

MS | A„, P», . 0*)+ A», ! ( 2 2 s ) P„-2„~i 0*) + 

= Mf*- 1 *!"!^,^)?.-*,^,^)*--"; (7) 
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in which the upper limit of p is governed by the same considerations as in 
the expression (5). 

6. Before formally enunciating the theorem that has been reached it may 
be verified, in its two cases, by applying y 2 to the expressions for the 
potential. 

Case 1.— r < /. 

Since 

V 2 <> 2 -~ 2 ?V cos ©'-fr' 2 )*™ 1 ' 2 = 2s(2s-l)(r a -2rr / cos «-fV 2 )*-^ 
we have 



t2 



/(/, 0') (r 2 -2rr cos B + r*Y~ V 2 A// d0' dr' 



o 

r27r 



= 2s(2s-l) /(/ J ^)(^-2r/cos@ + r ,2 ) s "' 3 / 2 ^ / ^ / ^/, (8) 

.J Jo 

and if (5) is the correct value of the potential the expression (8) may be 
deduced from (5) by writing therein s— 1 for s and multiplying the result 
by 2s(2s-l). 

The value of (8), thus obtained, is 

n = oo p /9 «_ 0\ 

2s(2s~l)M 2 JU» U w «a f ,,,-. 1 ,- 1 P M , p (^)r 11 ; (9) 

% = o p = o \ ^r I 

moreover, since 

y 2 P ms r n = 2s(2s-l)P n - aM - 1 r*- 8 l 

the operation of v 2 upon (5) yields 

M 2 S(o OT )a»-^,*- P 2^(2jp--l)P n -.2,jp-i(^)^*"" 2 . 
This may be written 

2s(2s—l)M ^ 2 o o)^-2^ S ^Pn-»2^-l(M)^ M ~ 2 , 

and, writing herein %+ 2 for n and _p + 1 for p, this is 

w = oo jo /2 s— -2\ 

2s(2s-l)M 2 2 '«„ U w ^2p, s ™^lPm^(^)^ 

ra = p ~ \ ^r ! 

agreeing with (9). 

7. Case2.— ?•>/. 

Similarly operating with y 2 upon (3) we obtain (8) as before ; and if (7) is 
the correct expression of the potential the value of (8) may be deduced 
from (7) by writing therein s— 1 for s, and subsequently multiplying by 
28(26*-!). 
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2s(2s-l)M>*-* 2 t ( g^ )A»„P^^i(/*)r-»; (10) 



The value of (8), thus obtained, is 

w = < » | ,2s— 2 

moreover since 

the operation of v 2 upon (7) yields 

n — co p /9 <? \ 

Mr3s - 3 2 s / » )A n , p (2s-2p)(2s-2p-l)V n - 2p ,,- p - 1 ( H ,)r-» 

n=0p=0 V'Jtfl 

which is 

n — co p /O q 0\ 

2s(2*-l)Mr»-» 2 2 ( 2 V A^pP.-^,-,-!^)*— 

agreeing with (10). 

8. The theorem being thus established, we may enunciate it in the 
following manner : — 

" Theorem. — The attracting body, of mass M and of uniform density, being 
a hollow shell whose exterior and interior bounding surfaces are surfaces of 
revolution about the axis of z; let the origin be taken within the interior 
bounding surface ; let 

be the distances from any points of the body to the points whose co-ordinates 
are z, 0, r, 6 respectively. 

Let SR ?) dM. = M£a n> s z n t 

or =M.%A ms z~^~ 2s+l \ 

the series which is convergent for the point z,Q being chosen. Then the 
potential at the point r, 9 is, when the series proceeding by positive powers of 
r is convergent, represented by 

2s- 1 
= MS < a n>s V ni + f o ) a n-2, s~l Pn> 1 + . . . + ( o rQ J a n~2p, s-p^mp + . . . > T n ) 

that is to say, it is obtained by substituting in the potential expression 

JxJLJitCtfi, gZ , 

r for z t and 

tor d ny s> 

the upper limits of p in the series written being reached when n— 2p is 



398 Major P. A. MacMahon and Mr. H. B. C. Darling. 

zero or unity, subject to the series not having already terminated by the 

vanishing of 

'2s \ 

Also when the series proceeding by negative powers of r is convergent, the 
potential is represented by 

2jR r q ctNL 



JxLZr < A», o JT n , § + f n ) An, i X n~*2j *~ 1 + . . . 



\ 2 7) I m ? *n-2p, s-p +*.. r T \ 



that is to say, it is obtained by substituting in the potential expression 
r for z, and 

ior a»> s 

the upper limit of p being determined as in the previous case." 

9. It will be noticed that in the series 

a ms Pn> O+f 2 ) a n-2,s-\Pntl J t «*»> 

A n ,o x n , s -{-( n ) A tt , 1 In— 2>*-i "I- •••> 

there is a remarkable correspondence. 
If the former be denoted by 

X (*> ?) 
the latter will be denoted by 

X ( p > A )- 

In the former series, every P has the first suffix n, 

10. We propose now to give the detailed results for the case of a ring of 
attracting matter, of radius c and of indefinitely small cross-section. 

For a ring of radius c, the potential at a point on the axis of z is 

M(c 2 +* 2 ;r 1/2 

and this may be expanded so that the numerical coefficients are functions 
P M , 5 with zero arguments. For convenience write 

x n, g \\j) == x ni s . 
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Then if z is less than c the expanded form is 



2 z i "1 /»\2n 



Hence a 2n , s - p 3nj / c -(a.-a.+i ) 

and the formula gives for the potential at the point r, 

Me*" 1 



Po> *' Po> o + ^ Ps, / Ps, o + ( o ) Po> s'-i P2, 1 f 72 



2s\ -r. / -^ . /2s\ -r, /^ !■ r 



+ *i ?4, s ?4j o + ( 2 ) ^ 2 > s ~ 1 -^ 4 ' * ~*~ ( 4 ) ^ >0 ' s ~ 2 ' ^ 4 > 2 



4 

4 



+ *j Pe> 5 I J 6j o + ( 2 ) P^> s-l' Pe? 1 + ( 4 ) P2, s-2 P6>2 + f *g j Po> s-3' p6, 3 f —$ 



& 



+ 



— M.C s Z < r2n>s P2n»0 + ( 9 ) P2w-2,s-l P2w>l + 



+ ( 2n) I ' () > s - n '^ 2n>n j [cj - 



Also if # be greater than c, the expanded form is 

m^sps,,/^ == m. 2 - 1 { p o,;+P2,/^) 2 +P4,/(?) 4 + 

Hence A 2n , , — P2* /c 2n 

and the formula gives for the potential at the point r, 6 

Mr 2 *' 1 



Po>o Poj s+ -{ ?2, o'P2, s+ ( 9 JP2, i/Pois-i r -3 



r 



^ S 1P. .'P. -_L/^ 5 W->. _ /- P. • I C 



+ < P4, o' P4j s + I 2 ) P 4 ' /P2j 5-1 + ( 4 ) P4, 2 Po, s- 



+ i P6>0 / P6,5+ ( 2 jP^j/P^s-l-f ( ^ )P6 J 2 / P2,s-2+(^ jPe, 3'P(),s-3 



C 6 
^6 



+ 



It will be noticed that the latter is obtainable from the former by inter- 
changing accented aud unaccented P functions, and also r and c. 

It will be clear that in the theory of attraction above considered the 
density need not be uniform throughout the mass. It will suffice if it is a 
function of r\ 6 r , and not of <£>'. 



vol. xcv. — a 2 K 



